Experimental Results from Heavy Ion Collisions at LHC 
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We review a subset of experimental results from the heavy-ion collisions at the Large Hadron 
Collider (LHC) facility at CERN. Excellent consistency is observed for various measurements like 
charged particle multiplicity density, second order azimuthal anisotropy and nuclear modification 
factor of charged hadrons across all the experiments at the LHC. Comparison to corresponding re- 
sults from the collisions at lower energy (centre of mass energy, ^snn = 200 GeV) at the Relativistic 
Heavy Ion Collider (RHIC) facility shows that system formed at LHC (^snn ~ 2.76 TeV) produces 
a larger number of particles per partcipating nucleon pair, attains higher energy density, has a larger 
system size, and lives for a longer time. Higher collision energy at LHC allows measurements to 
uncover interesting observations at kinematic regions in pseudorapidity and transverse momentum 
not previously measured in the heavy-ion collisions. Finally, we compare the LHC measurements to 
some model calculations to get a glimpse of the physical insight provided by these measurements. 

PACS numbers: 25.75.-q, 25.75.Ag, 25.75.Bh 
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I. INTRODUCTION 



The main goal of the high energy heavy-ion colli- 
sions is to study the phase structure of the Quantum 
Chromodynamic (QCD) phase diagram[l|-y]. One of 
the most interesting aspects of these collisions is the 
possibility of forming a phase of de-confined quarks 
and gluons, a system that is believed to have ex- 
isted in a few microsecond old Universe. First prin- 
ciple QCD calculations suggest that it is possible to 
have such a state of matter if the temperatures at- 
tained can be of the order of the QCD scale (^ 200 
MeV) [34a]- I^ laboratory, such temperatures could 
be attained by colliding heavy-ions at relativistic en- 
ergies. Further, at very high energies as encountered 
for heavy- ion collisions at LHC and RHIC, the de- 
confined phase could stay for long enough durations 
enabling to study the properties of the fundamental 
constituent (quarks and gluons) of any visible mat- 
ter. 

The results from heavy-ion collisions at RHIC 
have clearly demonstrated the formation of a de- 
confined system of quarks and gluons in Au-|-Au 
collisions at ^iii = 200 GeV [Mli- The pro- 
duced system exhibits copious production of strange 
hadrons, shovifs substantial collectivity developed in 
the partonic phase, exhibits suppression in high 
transverse momentum (px) hadron production rel- 
ative to p+p collisions and small fluidity as reflected 
by a small value of viscosity to entropy density ratio 
(rj/s). A factor of about 14 increase in energy for 
Pb+Pb collisions at LHC is expected to unravel the 
temperature dependence of various observables as 
well to provide an extended kinematic region in ra- 
pidity and pt where measurements in heavy-ion col- 



lisions did not exist previously. On the other hand, 
at RHIC a beam energy scan program is being un- 
dertaken to explore the other details of the QCD 
phase diagram |12| . 

In this review paper, we discuss a subset of results 
that have come out from LHC Pb+Pb collisions at 
^snn = 2.76 TeV. We have divided the discussion 
into three sections. In the second section we dis- 
cuss the consistency of various measurements among 
the three LHC experiments that have heavy-ion pro- 
grams: ALICE, CMS and ATLAS. In the third sec- 
tion, we make a comparative study between similar 
observables measured at lower energy collisions at 
RHIC to those from LHC. In doing this, we high- 
light the additional information that heavy-ion colli- 
sions at LHC bring compared to RHIC. In the fourth 
section, we present a comparison of various model 
calculations to the corresponding measurements at 
LHC. Finally we summarize our observations in the 
last section of the article. 



II. CONSISTENCY OF RESULTS AMONG 
LHC EXPERIMENTS 

A. Charged particle multiplicity 

One of the flrst measurements to come out of 
the heavy-ion collision program at LHC is the 
charged particle multiplicity per unit pseudorapidity 
in Pb+Pb colhsions at ^/^ = 2.76 TeV. Figure |I] 
shows the centrality (reflected by the number of par- 
ticipating nucleons, i^art, obtained from a Glauber 
model calculation [I3) dependence of dNch/drj at 
midrapidity for Pb+Pb collisions at ^/snn = 2.76 



10 







i Tf T i^ ■ 



• ALICE 2.76 TeV 
■ ATLAS 2.76 TeV 
A CMS 2.76 TeV 



% 50 100 150 200 250 300 350 400 

(N ) 

^ part ' 

FIG. 1: (Color online) Average charged particle mul- 
tiplicity density per unit pseudorapidity (dNch/drj) at 
midrapidity per participating nucleon (A'^part) pair plot- 
ted as a function of Apart for Pb+Pb collisions at ^snn 
= 2.76 TeV. The measurements are shown from AL- 
ICE [III, CMS [H and ATLAS ^ experiments. 



TeV from ALICE pj], CMS [T5| and ATLAS [Tg 
experiments. The error bars reflect statistical un- 
certainties. 

The ATLAS measurements oidNch/dri\,f=o are ob- 
tained over I77I < 0.5 using a minimum bias trigger 
with a central solenoid magnet off data set. The 



charged particles are reconstructed using two dif- 
ferent algorithms using the information from pixel 
detectors covering JTyj < 2.0. The A^part values are 
obtained by comparing the summed transverse en- 
ergy in the forward calorimeter over a pseudorapid- 
ity range 3.2 < \r]\ < 4.9 to a Glauber model simu- 
lation. The CMS results dNch/dri\rj=o are from the 
barrel section of the pixel tracker covering \ri\ < 2.5. 
The minimum bias trigger data set was in the mag- 
netic field off configuration so as to improve the ac- 
ceptance of low pt particles. The centrality determi- 
nation as in the case of ATLAS experiment are done 
using information from hadron forward calorimeter 
(2.9 < \r]\ < 5.2) and Glauber model simulations. 
The ALICE measurement uses a minimum bias data 
set from the silicon pixel detector (|7/| < 2.0). The 
centrality selection is carried out using signals from 
VZERO detectors (2 arrays of 32 scintillators tiles) 
covering the region 2.8 < 77 < 5.1 and -3.7 < 77 < - 
1.7, along with the corresponding Glauber modeling 
of the data. 



In-spite of difference in operating conditions and 
measurement techniques, the dNch/drj versus A'^part 
results for Pb+Pb collisions at ^snn = 2.76 TeV 
shows a remarkable consistency across the three ex- 
periments. 

B. Azimuthal anisotropy 
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Figure [2] shows the azimuthal anisotropy of produced charged particles {v„ = (cos(n(0 — ^n)))) as a 
function of pT for 30-40% Pb+Pb collisions at ^snn = 2.76 TeV. Here (f> is the azimuthal angle of the 
produced particles and VPn is the n*'* order reaction plane angle measured in the experiments. 

In the CMS experiment |17| the V2 measurements use the information from the silicon tracker in the region 
I77I < 2.5 with a track momentum resolution of 1% at pt ~ 100 GcV/c and kept within a magnetic field 
of 3.8 Tesla. The event plane angle (^2) is obtained using the information on the energy deposited in the 
hadron forward calorimeter. A minimum 77 gap of 3 units is kept between the particles used for obtaining 
\E'2 and V2- This ensures suppression of non-flow correlations which could arise for example from dijets. The 
event plane resolution obtained using three sub event technique varies from 0.55 to 0.84, depending on the 
collision centrality. The ATLAS experiment [l^ measured f„ using the inner detectors in the \ri\ < 2.5, 
kept inside a 2 Tesla fleld of superconducting solenoid magnet. The event planes are obtained using forward 
calorimeter information, with resolution varying from 0.2 to 0.85 depending on collision centrality. The 
ALICE experiment [l9| measured w„ using charged tracks reconstructed from the Time Projection Chamber 
(|7;| < 0.8), the event plane was obtained using information from VEZRO detectors kept at a large rapidity 
gap from the TPC. The momentum resolution of the tracks are better than 5%. 

A very nice agreement for V2, v^ and v^ values versus pt to a level of within 10% for most of the px range 
presented is found between all the experiments. 



C. Nuclear modification factor 

Figure [3] shows the nuclear modification factor for 
inclusive charged hadrons measured at midrapidity 
in LHC experiments for Pb+Pb collisions at ^snn 
= 2.76 TeV. The nuclear modification factor is de- 



fined as Raa = ^ — f''^ 1,1 ^iT — , here the overlap 

integral Tab = Nbinary/(T^neiasUc ■^ith Nbinary being 
the number of binary collisions commonly estimated 
from Glauber model calculation and daNN /dijd'^pT 
is the cross section of charged hadron production in 
p+p coUisions at ^^ =2.76 TeV. The result shows 
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FIG. 2: (Color online) v-n versus pr at midrapidity for 30-40% Pb+Pb collisions at ^snn = 2.76 TeV. The results 
are shown from different LHC experiments: CMS [13, ATLAS [Til and ALICE [li|. The bottom panels shows the 
ratio of the experimental data to a polynomial fit to the ALICE data. 
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FIG. 3: (Color online) Nuclear modification factor 
Raa of charged hadrons measured by ALICE [2y| and 
CMS [2l[ experiments at midrapidity for 0-5% most cen- 
tral Pb-Pb collisions at y'sNN = 2.76 TeV. The boxes 
around the data denote pT-dependent systematic uncer- 
tainties. The systematic uncertainties on the normaliza- 
tion are shown as boxes at -Raa = 1- 



that the charged particle production at high px in 
LHC is suppressed in heavy-ion coUisions relative to 
nuclcon-nucleon colhsions. 

The ALICE experiment [20] uses the Inner Track- 
ing System (ITS) and the Time Projection Chamber 
(TPC) for vertex finding and tracking in a minimum 
bias data set. The CMS experiment [2l| reconstructs 
charged particles based on hits in the sihcon pixel 
and strip detectors. In order to extend the statistical 
reach of the pt spectra in the highly pre-scaled min- 
imum bias data recorded in 2011 it uses unprescaled 
single-jet triggers. Both experiments take the value 

0f^fneZa.t.c=64±5mb. 

An excellent agreement for -Raa versus pt for 
charged hadrons in 0-5% central Pb-|-Pb collisions 



at .y/sNN ~ 2.76 TeV is observed between the two 
experiments. 

Having discussed the consistency of these first 
measurements in Pb-|-Pb collisions among different 
experiments, the major detectors used, acceptances 
and ways to determine centrality and event plane, we 
now discuss the comparison between measurements 
at RHIC and LHC heavy-ion collisions. 



III. COMPARISON OF LHC AND RHIC 
RESULTS 

In the first subsection, we discuss the energy de- 
pendence of basic measurements made in heavy-ion 
collisions. These include dNch/drj, {ttit) {fnT = 
\/Pt + mass^), Bjorken energy density (esj), hfe 
time of the hadronic phase (r/), system volume at 
the freczc-out, kinetic and chemical freeze-out con- 
ditions and finally the fiuctuations in net-charge dis- 
tributions. In the next subsection, we discuss about 
the energy dependence of pT integrated V2 which is 
a measure of collectivity, followed by comparisons 
of Vn versus px and flow fluctuations at RHIC and 
LHC. In the final subsection, we compare the nuclear 
modification factor for hadrons produced in heavy- 
ion collisions at RHIC and LHC. 



A. Bulk properties at freeze-out 

1. Multiplicity 

Figure [J] (top panel) shows the charged parti- 
cle multiplicity density at midrapidity [dNch/di]) 
per participating nuclcon pair produced in central 
heavy-ion collisions versus ^snn- We observe that 
the charged particle production increases by a fac- 
tor 2 as the energy increases from RHIC to LHC. 
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FIG. 4: (Color online) Top panel: dNch/drj per partci- 
pating nucleon pair at midrapidity in central heavy ion 
collisions as a function of ^snn- Bottom panel: Com- 
parison of dNch/drj per participating nucleon at midra- 
pidity in central heavy-ion collisions l isl. [iB . l22 - [3Cll | to 
correspondin g re sults from p+p{p) [3ll-[39l] and p{d)+A 
collisions [H.liO, EH . 



2. Average transverse mass and Bjorken energy 
density 
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FIG. 5: (Color online) Top panel: {mr) for charged 
pions in central heavy- ion collisions at midrapidity for 
AGS [13, m, SPS Ej,!!!!, RHIC ^\^ and LHC g^] 
energies. The errors shown are the quadrature sum of 
statistical and systematic uncertainties. Bottom panel: 
The product of Bjorken energy density, esj [4g], and 
the formation time (r) in central heavy ion collisions at 
mid-rapidity as a function of ^snn [49l - [54l ]. 




The energy dependence seems to rule out a logarith- 
mic dependence of particle production with ^snn 
and supports a power law type of dependence on 
^/suN. Figure |4] (bottom panel) shows the excess 
of dNch/dri/{Npart) in A-l-A colHsions [If 
l30l | over corresponding yields in p+p{p) p31l - 
p{d)+A collisions [23. l40l. l4l|. This observation also 
seen at RHIC persists at LHC but is proportionately 
larger at the higher energy collisions at the LHC. A 
power law fit to the p+p collision charged particle 
multiplicity density leads to a dependence ~ s"-^^, 
while those for A-l-A collisions goes as ^ s"-^^. This 
is a clear indication that A-f A collisions at RHIC 
and LHC are not a simple superposition of several 
p+p collisions, whereas the p-j-A collisions seem to 
closely follow this superposition principle. 



Figure[5]top panel shows the (tot) values for pions 
in central heavy-ion collisions as a function of y^snn- 
The (rriT) value increases with ^snn at lower AGS 
energies [42, 14311 ■ stays independent of ^snn for the 
SPS energies [4J, \^ and then tends to rise further 
with increasing ^snn at the higher beam energies 
of LHC. About 25% increase in {rriT) is observed 
from RHIC 



471. For a thermo- 



[33|, [461 to LHC 
dynamic system, (tot) can be an approximate rep- 
resentation of the temperature of the system, and 
dN/dy oc /^(^snn) rnay represent its entropy J55j . 
In such a scenario, the observations could reflect the 
characteristic signature of a phase transition, as pro- 
posed by Van Hove [Sg. Then, the constant value 
of (rriT) vs. y'sNN has one possible interpretation 
in terms of formation of a mixed phase of a QGP 
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FIG. 6: (Color online) Top panel: Product of the three 
pion HBT radii at kr (average transverse momenta of 
two pions) — 0.3 GeV/c for central heavy- ion collisions 
at AGS [sll, SPS m,!!!!, RHIC [13,111 and LHC [13 
energies. Bottom panel: The decoupling time extracted 
from Riongikr) for central heavy- ion collisions at midra- 
pidity at AGS, SPS, RHIC and LHC energies as a func- 
tion of (diVeh/dT))^''^ 



3. Volume and decoupling time 

The top panel of Fig. [5] shows the energy depen- 
dence of the product of the three radii (i?out, ^side 
and i?iong) obtained from pion HBT analysis. Here 
the "out" corresponds to the axis pointing along 
the pair transverse momentum, the "side" to the 
axis perpendicular to it in the transverse plane, and 
the "long" axis being along the beam (Bertsch-Pratt 
convention [6J,[63]). The product of the radii is con- 
nected to the volume of the homogeneity region at 
the last interaction. The product of the three radii 
shows a linear dependence on the charged-particle 
pseudorapidity density. The data indicates that the 
volume of homogeneity region is two times larger at 
the LHC than at RHIC. 

Further, within a hydrodynamic picture, the de- 
coupling time for hadrons (r^) at midrapidity can 
be estimated from the magnitude of radii i?iong- 
With i?2^„g = T]TK2{mT/T)/mTKi{mT/T), with 

niT = ■\/m^-|-fc|^, where m^ is the mass of the pion, 
T is the kinetic freeze-out temperature (taken as 120 
MeV) and Ki and K2 arc the integer order modified 
Bessel functions. The extracted r/ values for cen- 
tral heav y-io n collisions at midrapidity at AGS 57| , 
SPS [H 111, RHIC [63,[6l| and LHC ^ energies 
are shown as a function of cube-root of dNch/drj in 
the bottom panel of Fig. El We observe that tj 
scales linearly with {dNch/dviY/'^ and is about 10 
fm/c at LHC energies. Note the value is about 40% 
larger than at RHIC. However there are few caveats 
which needs to be kept in mind, the above expres- 
sion ignores transverse expansion of the system, fi- 
nite chemical potential for pions and uncertainties 
associated with freeze-out temperature. These ef- 
fects are to some extent reflected in the uncertainties 
shown on the data points. 



and hadrons during the evolution of the heavy-ion 
system. The energy domains accessed at RHIC and 
LHC will then correspond to partonic phase while 
those at AGS would reflect hadronic phase. How- 
ever, there could be several other effects to which 
(tot) is sensitive, which also need to be understood 
for proper interpretation of the data [55| . 

Figure[5]bottom panel shows the product of the es- 
timated Bjorken energy density [esj = A^dEx/dy; 
A± [45 is the transverse overlap area of the nuclei 
and Et is the transverse energy) and formation time 
(r) as a function of -^snn |49l454l |. The energy den- 
sity at LHC seems to be a factor 3 larger compared 
to those attained at RHIC. 



4- Tem,perature and collective flow velocity 

The hadron yields and spectra reflects the proper- 
ties of the bulk matter at chemical and kinetic freeze- 
out, respectively. Generally, the point at which 
the inelastic collisions ceases is called the chemical 
freeze-out and the point where even the elastic col- 
lisions stop is called the kinetic freeze-out. 

The transverse momentum distribution of differ- 
ent particles contains two components, one random 
and the other collective. The random component 
can be identified with the temperature of the system 
at kinetic freeze-out (Tkin). The collective compo- 
nent, which could arise from the matter density gra- 
dient from the center to the boundary of the fireball 
created in high energy nuclear collisions is called col- 
lective flow in transverse direction ((/3)). Using the 
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FIG. 7: (Color online) Kinetic freeze-out temperature 
and collective flow velocity in central heavy-ion collisions 
as a function of collision energy. [33 . |47| . l65l - [6a | 



assumption that the system attains thermal equi- 
librium, the blast wave formulation can be used to 
extract Tkin and (/3) . These two quantities are shown 
in Fig. [7] versus ^snn |3j, |43, [65| - |68| . For beam en- 
ergies at AGS and above, one observes a decrease 
in Tkin with y^snn- This indicates that higher the 
beam energy, longer is the interactions among the 
constituents of the expanding system and lower the 
temperature. From RHIC top energy to LHC there 
seems to be however a saturation in the value of 
Titin- In contrast to the temperature the collective 
flow increases with increase in beam energy rapidly, 
reaching a value close to 0.6 times the speed of light 
at the LHC energy. 

Figure [5] shows the chemical freeze-out tempera- 
ture (Teh) versus the baryon chemical potential (/ie) 
in central heavy- ion collisions [3j, |43, [69l - l76j . These 
quantities are obtained by fitting the particle yields 
to a statistical model assuming thermal equilibrium 
within the framework of a Grand Canonical ensem- 
ble. There are two values of temperature quoted 
for LHC energies. A Teh value of about 164 MeV 
and fixed /is value of 1 MeV seems to reproduce 
the multi-strange ratios (involving S and fl) quite 
well, but were observed to miss the data for p/ir 
and A/tt. On the other hand, the statistical ther- 
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FIG. 8: (Color online) Chemical freeze-out tempera- 
ture versus baryon chemical potential in central heavy- 
ion collisions [33. 1471. [69l - [76l [ . The curve corresponds to 
model calculations from Ref. [69l.l7(]|. 



mal model prediction with Teh — 152 MeV and fixed 
/iB = 1 MeV fits the measured p/ir and A/tt ratios 
better but misses the ratios involving multi-strange 
hadrons [73] ■ This issue is not yet resolved, with 
being possibly related to hadronic final state inter- 
actions [78| . The curve corresponds to model calcu- 
lations from Ref. ^^. 



5. Fluctuations 
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FIG. 9: (Color online) Energy dependence of net-charge 
fluctuations about midrapidity in central heavy-ion colli- 
sions at SPS [zil, RHIC iOl and LHC [H energies. Also 
shown are the expectations from a hadron resonance gas 
model and for a simple QGP picture [82 |. 



The fluctuation in nct-chargc depends on the 
square of the charge state of the system. They are 
good indicators of phase transition, where the par- 
tonic phase have constituents with fractional charges 
while the hadronic phase has constituents with in- 
tegral units. Hence, net-charge fluctuations are ex- 
pected to be smaller if the system underwent a phase 
transition and freeze-out is close to the phase tran- 
sition temperature so that the signals survive. An 
experimental measure of net-charge fluctuations is 
defined as v{+-,dyn) - <^+^^+-^» ^ {n.^n.-i)) 
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- 2 /jv )(N ) ■• ■^here {N-) and (A^+) are average neg- 
ative and pos itive charged particle multiplicity, re- 
spectively |83j . 

Figure [S] shows the product of v{-\ — ,dyn) and 
{Nch) (average number of charged particles) as a 
function of ^snn |79l - [8lj . We find that this fluc- 
tuation observable rapidly decreases with ^/stssn and 
approaches expectation from this observable for a 
simple QGP-like scenario [S^] as we move from RHIC 
to LHC energies. Given that several other observ- 
ables already indicate that a hot and dense medium 
of color charges has been formed at RHIC and 
LHC, this result may indicate that the observable 
i'{-\ — ,dyn) is not sensitive enough to QGP physics 
or the process of hadronization washes out the QGP 
signal for this observable. It may be also noted that 
the models results does not incorporate the accep- 
tance effects and docs not consider any dynamic evo- 
lution of the system like for example the dilution of 
the signals in the hadronization process. 
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FIG. 10: (Color online) Transverse momentum inte- 
grated V2 close to midrapidity for charged (Z = 1) parti- 
cles for collision centralities around 20-30% as a function 
of center of mass energy. 



tion of beam energy is observed. 



The experimental 
data used are from FOPI [85|,j8dl, EOS, E895 [13, 
E877 pi, CERES [8ijNA49 (qJ, STAR [lH, PHO- 
BOS [al, PHENIX|3, ALICE [H, ATLAS [H 
and CMS [131 experiments. Charged particles are 
used for LHC, RHIC, CERES and E877 experi- 
ments, pion data is used from NA49 experiment, 
protons results are from EOS and E895 experiment 
and FOPI results are for all particles with Z=l. 



B. Azimuthal anisotropy 



2. Azimuthal anisotropy coefficients versus transverse 
momentum 



1. Energy dependence of Pt integrated V2 

Figure [TU] shows the pt integrated V2 close to 
midrapidity for charged particles for collision cen- 
tralities around 20-30% as a function of center of 
mass energy. We observe that there is an increase 
in magnitude of V2 by about 30% from top RHIC 
energy {^/suN = 200 GeV) to LHC energy {^/sm^ = 
2.76 TeV). This needs to be viewed within the con- 
text of a similar magnitude of increase in (tot) of 
pions from RHIC to LHC energies. The increase of 
V2 beyond beam energy of 10 AGeV is logarithmic 
in ^snn- This is expected to be determined by the 
pressure gradient-driven expansion of the almond- 
shape fireball produced in the initial stages of a 
non-central heavy-ion collision [8J]. While for V2 
measured at lower beam energies, the dependences 
observed is due to interplay of passing time of spec- 
tators and time scale of expansion of the system. 
A preference for an in-plane emission versus out-of- 
plane ( "squeeze-out" ) pattern of particles as a func- 



Figure ITlT top) shows the comparison of V2{pt), 
V3{pt) and V4{pt) at RHIC^HENIX experi- 
ment [H) and LHC (ALICE [9^) at midrapidity 
for 30-40% Au-fAu and Pb-hPb collisions respec- 
tively. The bottom panel of this figure shows the 
ratio of LHC and RHIC results to a polynomial fit 
to the LHC data. The w„(pt) measurement tech- 
niques are similar at RHIC and LHC energies. One 
observes that at lower px { < 2 GeV/c) the V2{pt) 
and vz{pt) are about 10-20% smaller at RHIC com- 
pared to the corresponding LHC results. However at 
higher pT the results arc quite similar. The W4(pt) 
seems higher at RHIC compared to LHC. 

One of the most striking observations to come out 
from RHIC is the number of constituent quark [uq) 
scaling of V2{pt) for identified hadrons. The ba- 
sis of such a scaling is the splitting of V2(j>t) be- 
tween baryons and mesons at intermediate pt (2-6 
GeV/c). This is shown in the bottom panels of the 
bottom Fig. 111! Such a splitting between baryon and 
meson V2{pt) is also observed at intermediate pr at 
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FIG. 11: (Color online) Top: Comparison ol v„{pt) 
at midrapidity for 30-40% collision centrality at RHIC 
(Au-HAu collisions at ^^^s^ = 200 GeV from PHENIX 
experiment [93] ) and at LHC (Pb-fPb collisions at -^snn 
= 2.76 TeV from ALICE experiment [qB]). The bottom 
panel shows the ratio of Vn at LHC and RHIC. Bottom: 
V2 versus pt and vijuq versus pr /riq for pion and pro- 
tons at midrapidity for 10-20% collision centrality from 
Au-I-Au collisions at y/a^ = 200 GeV (PHENIX exper- 
iment [H) and Pb-I-Pb collisions at y/s^ = 2.76 TeV 
(ALICE experiment [H). 



LHC energies (seen in the top panels of the bottom 
Fig. [TT|) . However the degree to which Uq scaling 
holds could be different at RHIC [H and LHC [if 
energies. The Uq scaling is much more closely fol- 
lowed at RHIC compared to LHC. It may be noted 
that there are several factors which could dilute such 
scalings, which includes energy dependence of radial 
flow, an admixture of higher Fock states and con- 
sideration of a realistic m omen tum distribution of 
quarks inside a hadron [93, llOOf . The observation of 
the baryon-mcson splitting is commonly interpreted 
as due to substantial amount of collectivity being 
generated in the dc-confincd phase. Another impor- 
tant feature is that at both RHIC and LHC energies 
a clear hydrodynamic feature of mass dependence of 
V2[pt) is observed at low pT {< 2 GeV/c). 

Figure [T^ shows the charged hadron V2{pt) for 
30-40% collision centrality in Au+Au collisions at 
= 200 GeV and Pb+Pb collisions at -v/snn = 



2.76 TeV for \r]\ < 1 [17|. This figure demonstrates 
the kinematic reach for higher energy collisions at 
LHC relative to RHIC. LHC data allows us to study 




CMS Pb-Pb ^s^;^ = 2.76 TeV 
STAR Au-Au f^ = 200 GeV 
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FIG. 12: (Color online) Comparison of V2{pt) at midra- 
pidity for 30-40% collision centrality at RHIC (Au+Au 
collisions at ^snn = 200 GeV from STAR experiment) 
and at LHC (Pb-|-Pb collisions at y/s^ = 2.76 TeV 
from CMS experiment [l3l). The shaded band about 
CMS data point are systematic errors and vertical lines 
represent statistical errors. 



the V2{pt) in the px range never measured before 
in heavy-ion collisions. The V2{pt) ^ for pT > 
40 GeV/c suggesting that those particles must have 
been emitted very early in the interactions, when 
collective effects had not set in. These high trans- 
verse momentum data is useful to understand the 
effects of the initial geometry or path-length depen- 
dence of various properties associated with parton 
modification inside the hot QCD medium. In addi- 
tion, it also provides significantly improved precision 
measurement of V2 for 12 < pT < 20 GeV/c. 



3. Flow fluctuations 

Fluctuations in azimuthal anisotropy coefficient 
V2 have gained quite an attention in recent times. 
In particular, the measurement of event- by- event 
V2 fluctuations can pose new constraints on the 
models of the initial state of the collision and its 
subsequent hydrodynamic evolution. In extracting 
cvcnt-by-event V2 fluctuations one needs to sepa- 
rate non-flow effects and so far there is no direct 
method to decouple V2 fluctuations and non-flow 
effects in a model independent way from the ex- 
perimental measurements. However, several tech- 
niques exists where the non-flow effects can be mini- 
mized, for example, flow and non-flow contributions 
can be possibly separated to a great extent with 
a detailed study of two particle correlation func- 
tion in A(f> and its dependence on rj and A77. Here 
we discuss another technique to extract and com- 
pare the V2 fluctuations at RHIC and LHC. We as- 
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sunic that the difference between f2{2} (two par- 
ticle cumulant) and U2{4} (four particle cumulant) 
is dominated by V2 fluctuations and non-flow effect 
is negligible for W2{4}. Then the ratio Rvi2-4) ~ 
V(«2{2}2 _ i,2{4}2)/(t;2{2}2 -|- W2{4}2) cau be con- 
sidered as an upper limit for V2 fluctuations in the 
data. Figure [T3] shows the i?„(2-4) as a fun ction 
of collision centrality and dNch/drj for RHIC |lOl| 
and LHC [O^ energies. The centrality dependence 
of i?„(2-4) at RHIC or LHC as seen in Fig. [T3l could 
be an interplay of residual non-flow effects which in- 
creases for central collisions and multiplicity fluctua- 
tions which dominates smaller systems. It is striking 
to see that Ryi2-i) when presented as a function of 
% cross section is similar at RHIC and LHC suggest- 
ing it reflects features associated with initial state of 
the collisions for example the event- by- event fluctu- 
ations in the eccentricity of the system. But when 
presented as a function of dNch/dij it tends to sug- 
gest a different behaviour for central most collisions 
at RHIC. 



200 400 600 800 1000 1200 1400 1600 
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FIG. 13: (Color online) The ratio i?„(2-4) = 
^{v2{2}'2 - U2{4}2)/(«2{2}2 + U2{4}2), a mcasuie of 
upper limit of V2 fluctuations plotted as a function 
of collision centrality (top panel) and average number 
dNch/dr) (bottom panel) for RHIC (ST AR e xperiment: 
Au-I-Au collisions at ^/s^ = 200 GeV [loH ') and LHC 
(ALICE: Pb+Pb collisions at ^/si57 = 2.76 TeV ^) at 
midrapidity. The bands reflect the systematic errors. 



Nuclear modification factor 



J 



Figure [T3] shows the Raa of various particles produced in heavy- ion collisions at RHIC and LHC. In 
Fig. [Mia), we observe that the shape of the i?AA versus pT of charged hadrons at RHIC and LHC [20, l2l| 
are very similar for the common p^ range of measurements. The values i?AA at RHIC are higher compared 
to those at LHC energies up to pt < 8 GeV/c. The higher kinematic reach of LHC in pT allows us to see the 
full pt evolution of i?AA in high energy heavy-ion collisions. All these suggests the energy loss of partons in 
the medium formed in heavy-ion collisions at LHC energies is perhaps larger compared to those at RHI C. In 
Fig. rMT b). we observe that the nuclear m odifi cation factors for d-l-Au collisions at y^snn = 200 GeV jl02| 
and p-f Pb collisions at ^/snH = 5.02 TeV |l04l | are greater than unity for the pr > '2 GeV/c. The values for 
RHIC are slightly larger compared to those for LHC. A value greater than unity for t he nuclear modification 
factor in p{d)+A collisions are generally interpreted as due to Cronin effect (lOSlllOq . However several other 
physics effects could infiuence the magnitude of the nuclear modification factor in p{d)+A collisions, which 
includes nuclear shadowing and gluon saturation effects. But the fact that the nuclear modification factor 
in p{d)+A collisions are not below unity, experimentally strengthens the argument for formation of a hot 
and dense medium of color charges in A-l-A collisions at RHIC and LHC. In Fig. [Hfc). we show the i?AA 
of particles than do not undertake strong interactions and some of th em a re most like ly formed in the very 
early stages of the collisions. These particles (photon [l07lll08l |. W^ |l09f and Z J110l | bosons) have a i?AA 
~ 1, indicating that the i?AA < 1 observed for charged hadrons in A-f-A collisions are arising due to strong 
interactions in a dense medium consisting of color charges. 



10 




'l"'l' 



' I " ' I " ' I ' " I "i' 



I ' ' ' I ' ' ' I ' 

CMS PbPb (i^ = 2. 

• W 
AW* 
T W 
*Z° 



CMS PbPb (0-5%) ^=2.76 TeV 
ALICE PbPb(0-5%) ^=2.76TeV 
STAR AuAu(0-5%) ^=0.2TeV 
PHENIX AuAu (0-10%) ^= 0.2TeV 



.*++■ *H 




• ALICE pPb^ = 5.02TeV 
ir STAR dAu)(s;;„ = 0.2TeV 




I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I 



(C) 



ffl 



Isolated photon 

■ CMS PbPb (0-10%) 

PHENIX AuAu (0-5' 

1 ■ ■ ■ I ■ ■ ■ I 






76TeV 
0.2 TeV 



_L. 



2 3 456 10 20 30 



P, (GeV/c) 



100 2 4 6 8 10 12 14 16 18 20 22 

p^ (GeV/c) 



20 40 60 80 100 120 



P, (GeV/c) 



FIG. 14: (Color online) (a) Nuclear modification factor Raa of charged hadrons measured by ALICE [20| and 
CMS [2l| experiments at midrapidity for 0-5% most central Pb-fPb collisions at ^snn ~ 2.76 TeV. F or co mparison 
shown the -Raa of charged hadrons at midrapidity for 0-5% most central collisions measured by STAR [l02| and i?AA 
of tt" at midrapidity for 0-10% most central collisions measured by PHENIX [l03l ] for Au-|-Au collisions at -^/snn = 
200 GeV. The boxes around the data denote pT-dependent systematic uncertainties. The systematic uncertainties 
on the normalization are shown as boxes at -Raa ~ 1. (b) Comparison of nuclear modification factor for charged 
hadrons versus pr at midrapidit y for minimum bias collisions in d+Au collisions at ^snn ~ 200 GeV [l021 | and p-|-Pb 
collisions at ^snn ~ 5.02 TeV [lOll ]. (c) The nucle ar modification fact or versus pr for isolated photons i n ce ntral 
nucleus-nucleus coUisions at y/s^ = 200 GeV [TotI and 2.76 TeV [lo3 l. Also shown are the Raa of W* [l09l ] and 
Z bosons [110| | at LHC energies. 



IV. COMPARISON TO MODEL 
CALCULATIONS 

In this section, we compare some of the experi- 
mental observables discussed above with correspond- 
ing model calculations. This helps us to interpret 
the data at both RHIC and LHC energies. The ob- 
servables chosen are the charged particle multiplicity 
density, ratio of kaon to pion yields as a function of 
beam energy, and transverse momentum dependence 
of V2 and -Raa • 



A. Charged particle multiplicity density and 
particle ratio 

Figure [T51 compares the measured charged particle 
pseudorapidity density at RHIC (0.2 TeV) and LHC 
(2.76 TeV) energies to various model calculations. 

Empirical extrapolation fro m lo wer energy data 
(named "Busza" in the figure) |lll| significantly un- 
der predicts the measurement at LHC energies. A 
simple power-law growth of charged-particle mul- 
tiplicities near midrapidity in central Au-|-Au col- 
lisions seem to be followed up to RHI C en ergies 
(named as "Barshayet. al" , in the figure) jll2| . Per- 
turbativG QCD-inspired Monte Carlo event genera- 
tors, the HIJING model wi thou t jet quenching (113J |. 
the Dual Parton Model ^1^ (named "DPMJET 
III" in the figure), or the Ul tra- r elativistic Quantum 
Molecular Dynamics model jll5| (named "UrQMD" 



in the figure) are consistent with the measurement. 
The HIJING model results without jet quenching 
was also consistent with the RHIC measurements. 
The semi-microscopic models like LEXUS are suc- 
cessful in explaining the observed multiplici ty a t 
RHIC (named as "Jeon et. al" in the figure) J116l |. 
Models based on initial-state gluon density satura- 
tion have a range of prediction s depend ing on the 
specifics of the implementation |ll7l4l21 1 . The best 
agreement with LHC data happens for model as de- 
scribed in ( named as "Kharzeevl" and "Armesto" in 
the figure) [ll9l Il2ll | . Similar are the conclusions for 
RHIC energy from these models. The prediction of 
a hybrid model based on hydrodynamics and satu- 
ration of final-state phase space of scat tered partons 
(named as "Eskola" in the figure) (122| is slightly on 
higher side compared to the measurement at LHC. 
But such a mo del s eems to do a reasonable job for 
RHIC energies jl23j . Another hydrodynamic model 
in which multiplicity is scaled from p+p collisions 
over predicts t he m easurement (named as "Bozek" 
in the figure) |l24J . Models incorporating scaling 
based on Landau hydr odyn amics (named as "Sark- 
isyan" in the figiure) [l25l | and based on modified 
PYTHIA and hadroni c re-s cattering (named as "Hu- 
manic" in the figure) jl26| under predicts the mea- 
surement at LHC energy. At RHIC energies, models 
considering minijet production in ultra-relativistic 
heavy-ion collisions by taking semi-hard parton re- 
scatterings explicitly into account under predicts the 
multiplicities (named as "Accardi et. al in the fig- 
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FIG. 15: (Color online) Comparison of dNch/dr] mea- 
surement at midrapidity for central heavy-ion collisions 
at RHIC and LHC with model predictions. 



ure) [127l |. It is also seen at R HIC energies that 
models bas ed on string fusion |l28l |. dual string 
model [l29l | seem to work well, whereas those based 
on heavy ion cascade LUCIFER model |13C| under 
predicts the data. 

Figure [TBI shows the (rfiVch/(i?/)/((iVpart)/2) ver- 
sus (iVpart) for Pb-fPb colhsions at ^snn = 2.76 
TeV [iJl. Also shown are the corresponding RHIC 
results scaled up by a factor 2.15. Remarkable simi- 
larity is observed in the shape of the distributions 
at RHIC and LHC energies. Particle production 
based on saturation model explains the tr ends nicely 
(named as "Albacete et al" in the figure) [l3ll | . Sim- 
ple fit to the data using a power law form for the 
(A^part) also explains the measurements. In addition 
a functional form inspired by the detailed shape of 
pseudorapidity distribution of charged particle mul- 
tiplicity distributions at RHIC [331 explains the cen- 
trality trends nicely. 

Strangeness production in heavy-ion c ollisi ons is 
a classic signature for formation of QGP |l32| . The 
particle yield ratio K/tt could reflect the strangeness 



FIG. 16: (Color online) Centrality dependence of 
(diVch/d77)/((iVpart)/2) on (Np^rt) for Pb-(-Pb collisions 
at y^SNN = 2.76 TeV and Au+Au collisions at .^snn 
= 200 GeV. The RHIC results are scaled up by a fac- 
tor of 2.15. Also shown are comparisons to theoretical 
model calculations and some parametrization based on 
detail shape of dNch/drj distributions at RHIC [3y] and 

(iVpart). 
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FIG. 17: (Color online) Energy dependence of K /it 
ratio for central collisions at midrapidity. Errors are sta- 
tistical and systematic added in quadrature. Results 
are also compare d with various theoretical model pre- 
dictions [l3i4l3il. 



enhancement in heavy-ion collisions with respect to 
the elementary collisions. Figure [T7] shows the en- 
ergy dependence of K^/n^ ratio for central colli- 
sions at midrapidity. It will be interesting to see 
which model explains such an impressive collection 
of systematic data on K/n ratio. Figure [T7] also 
shows the energy dependence of K/tt ratio from var- 
ious theoretical model calculations. The energy de- 
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pcndcncc of /v +/7r+ ratio has been interpreted usin g 
the Statistical Model of Early Stage (SMES) pSJ . 
The model predicts first order phase transition and 
the existence of mixed phase around beam energy 
of 7-8 GeV. The SHM or Statistical Hadronization 
Model |l34l | assumes that the strong interactions sat- 
urate the particle production matrix elements. This 
means that the yield of particles are controlled pre- 
dominantly by the magnitude of the accessible phase 
space. The system is in chemical non-equilibrium 
for ^/SNN < 7.6 GeV, while for higher energies, 
the over-saturation of chemic al o ccupancies is ob- 
served. The Statistical Model jl35l | assumes that the 
ratio of entropy to T^ as a function of collision en- 
ergy increases for mesons and decreases for baryons. 
Thus, a rapid change is expected at the crossing 
of the two curves, as the hadronic gas undergoes 
a transition from a baryon-dominated to a meson- 
dominated gas. The transition point is character- 
ized by r=140 MeV, fj.B= 410 MeV, and ys]^=8.2 
GeV. In the Thermal Model |l36l | , the energy depen- 
dence of K /tt^ is studied by including cr- meson, 
which is neglected in most of the models, and many 
higher mass resonances (m > 2 GeV/c^) into the res- 
onance spectrum employed in the statistical model 
calculations. The hadronic non-equilibrium Kinetic 
model jl37l | assumes that the surplus of strange par- 
ticles is produced in secondary reactions of hadrons 
generated in nuclear collisions. Then the two impor- 
tant aspects are the available energy density and the 
lifetime of the fireball. It is suggested that these two 
aspects combine in such a way so as to show a sharp 
peak for the strangeness-to-entropy or K/tt ratio as a 
function of beam energy. In the Hadron Reso nanc e 
Gas and Hagedorn model (HRG+Hagedorn) jl38j |. 
all hadrons as given in PDG with masses up to 2 
GeV/c^ are included. The unknown hadron res- 
onances in this model are included through Hage- 
dorn's formula for the density of states. The model 
assumes that the strangeness in the baryon sector 
decays to strange baryons and does not contribute 
to the kaon production. The energy dependence 
of K/tt ratio seems to be best explained using 
HRG+Hagedorn model. 

This systematic measurement of K/tt ratio reveals 
two interesting information: (a) Around y/swN ~ 
8 GeV, the K^ /tt^ ratio shows a peak while the 
K^ /tt^ ratio increases monotonically indicating the 
role of the maximum baryon density at freeze-out 
around this collision energy and (b) for ^snn > 100 
GeV, pair production becomes the dominant mecha- 
nisim for K^ production, so both the ratios K^ /tt^ 
and K~ /tt~ approach the value of 0.16. Taking into 
account of different masses between pions and kaons, 
this asymptotic value corresponding to a tempera- 
ture of the order of 160 MeV. 



B. Azimuthal anisotropy 
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FIG. 18: (Color online) The azimuthal anisotropy pa- 
rameter V2 measured in non-central heavy-ion collisions 
at midrapidity for RHIC and LHC energies. For com- 
parison, shown are the various theoretical calculations 
based on hydrodynamic and transport approaches (see 
text for details). 



The azimuthal anisotropy parameter V2 measured 
at RHIG and LHG provides an unique opportunity 
to study the transport properties of the fundamen- 
tal constituents of any visible matter - a system of 
quarks and gluons. Further, it provides an opportu- 
nity to understand whether the underlying dynamics 
of the evolution of the system formed in the collisions 
is g overned by macroscopic hydrodynamics ll3SJ - 
Il4l| or by microscopic transport approach [142| . 
Figure [TH] shows the V2 versus pr for 30-40% col- 
lision centrality Au-|-Au and Pb-|-Pb collisions at 
midrapidity for ^/s^ = 200 GeV and 2.76 TeV, 
respectively. The measurements are compared to a 
set of model calculations based on hy drodynam ic ap- 
proach (including THERMINATOR [Hi [HI) and 
another set of calculations based on transport ap- 
proach. It is observed that hydrodynamic based 
models explain the V2 measurements both at RHIC 
and LHC energies. Transport based mode ls includ- 
ing partonic interactions (like AMPT jl42l | ) also ex- 
plain the V2 measurements. However, those trans- 
port models which do not incorp orate partonic in- 
teractions like UrQMD [lil Il46j fail to explain the 
data. The model comparison also reveals that data 
seem to suggest a high degree of fluidity reflected 
by a small value of shear viscosity to entropy den- 
sity ratio (77/s) < 0.2 being supported by the mea- 
surements. A more detail comparison of the model 
calculations with various order azimuthal anisotropy 
parameter w„ would in near future give us a more 
quantitative picture of the temperature (or energy) 
dependence of transport coefficients of the system 
formed in the heavy-ion collisions. 
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FIG. 19: (Color online) Measurements of the nuclear modification factor Ra a in central heavy-ion collisions at two 
different center-of-mass energies, as a function of pT, for pions (vr '°) [I47l . ll48 | and charged hadrons [20l . [2ll |. compared 
to several theoretical predictions (see text). The error bars on the points are the statistical uncertainties, and the 
boxes around the data points are the systematic uncertainties. Additional absolute normalization uncertainties of 
order 5% to 10% are not plotted. The bands for several of the theoretical calculations represent their uncertainties. 

The nuclear modification factor (Raa) is an observable used to study the structure of strongly interacting 
dense matter formed in heavy-ion collisions. Here we discuss the observation of Raa < 1 at high px seen at 
RHIC and LHC by comparing to models within perturbative QCD (pQCD) based formalisms. In this picture, 
the high pT hadrons are expected to originate from the fragmentation of hard partons (hard scattering scales 
large than QCD scales of 200 MeV). The hard partons lose energy through interactions with the hot and 
dense medium, which get reflected in the observed values of Raa- The processes by which they could loose 
energy includes radiative energy loss and elast ic energy loss. For a more elaborate discussion on these models 
we refer the reader to the review article jl49{ . 

In Fig. I19[ we show a comparison between experimentally measured Raa versus pt at LHC and RHIC 
energies to corresponding pQCD based model calculations. All theoretical formalisms requires a microscopic 
model of the medium to set the input parameters for the energy loss calculation. These parameters for 
example are denoted as (q), the transport coefficient of the medium or the gluon number density dN^ /dy 
per unit rapidity. The parameter Pesc on the other hand, reflects the strength of elastic energy loss put in 
the model calculations. Without going into deeper theoretical discussions of each mo del, we refer the readers 
to the following related pubheations: PQM ^150], GLV [lll|, ASW [HI, YaJEM [Hl, WHDG ^ and 
ZOWW jl55j |. However, for completeness and to elucidate the approach taken in the model calculations, we 
briefly mention as a example that two formalisms: the GLV approach named after their authors Gyulassy, 
Levai and Vitev and ASW approach named after the corresponding authors Armesto, Salgado and Wiede- 
mann model, the medium as separated heavy static scattering centers with color screened potentials. Where 
as in some other formalism a more precise definition of the medium is considered like being composed of 
quark gluon quasi-particles with dispersion relations and interactions given by the hard thermal loop effective 
theory. 

We observe that most models predict the pr dependence of Raa well for collisions both at RHIC and LHC 
energies. The models specially capture the generally rising behavior of Raa that is observed in the data 
at high Pt for the LHC energies. The magnitude of the predicted slope of Raa versus pT varies between 
models, depending on the assumptions for the jet-quenching mechanism. The models do not seem to indicate 
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the need for larger values of inedium density parameters in the calculation to explain the Raa at RHIC and 

LHC for the common kinematic range. However, it seems they may require a high medium density to explain 

the values of Raa for pr > 20 GeV/c. 

I 



V. SUMMARY 

In summary, the results on multiplicity density in 
pseudorapidity, HBT, azimuthal anisotropy, nuclear 
modification factor from LHC experiments indicate 
that the fireball produced in the nuclear collisions is 
hotter, lives longer and expands to a larger size at 
freeze-out compared to lower energies. These results 
also confirm the formation of a de-confined state of 
quarks and gluons at RHIC energies and provide 
a unique kinematic access in the measurements to 
study in detail, the properties (such as transport co- 
efficients) of this system of quarks and gluons. 

In this review, we first showed that the first set 
of measurements made by the three LHC experi- 
ments with a heavy-ion program, ALICE, CMS and 
ATLAS, show a high degree of consistency. These 
measurements includes, centrality dependence of 
charged particle multiplicity, azimuthal anisotropy 
and nuclear modification factor measurements ver- 
sus transverse momentum. Next we discussed the 
comparison of various measurements made at RHIC 
and LHC energies. LHC measurements of dNch/dr] 
clearly demonstrated the power law dependence of 
charged particle multiplicity on the beam energy. 
The LHC results also reconfirmed the observation at 
RHIC that particle production mechanism is not a 
simple superposition of several p+p collisions. The 
{rriT), EBj, freeze-out volume, decoupling time for 
hadrons, (W2) and (/3) are larger at LHC energies 
compared to RHIC. This is even though the freeze- 
out temperatures are comparable. The net-charge 
fluctuations measure rapidly approach from RHIC 
to LHC towards a simple model based calculation for 
QGP state. Although the senstivity of this observ- 



able for a heavy-ion system as well as lack of proper 
modeling of the heavy-ion system theoretically for 
such an observable needs careful consideration. The 
V2 fluctuations as a function of centrality fraction 
have a similar value at both RHIC and LHC. This 
reflects their sensitivity to initial state effects. Just 
like the d-l-Au and direct photon Raa measurements 
experimentally proved that the observed Raa < 1 
for charged hadrons is a final state effect at RHIC, so 
also the p-l-Pb, direct photon, W^ and Z° Raa mea- 
surements at LHC showed that the observed i?AA < 
1 is indeed due to formation of a dense medium of 
colored charges. All these conclusions were further 
validated by the comparison of several observables 
to corresponding model calculations in the previous 
section. This includes the finding that the fluid at 
LHC still shows a high degree of fluidity reflected by 
a small value of shear viscosity to entropy density 
ratio. 

M easureme nts related heavy quark produc- 
tion |l56l4l58i. dilep ton production, jet-hadron cor- 
relations Il59l . II6OI I and higher order azimuthal 



anisotropy 16ll . Il62| which are now coming out of 
both RHIC and LHC experiments will provide a 
much more detail characterization of the properties 
of the QCD matter formed in heavy-ion collisions. 
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